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1 Introduction 

The problem of locating sensors on the network to determine traffic flow volumes has been object in the 
past few years of growing interest. Problems in this class differentiate according to the type of sensors 
that need to be located (counting sensors, path-ID sensors, vehicle-ID sensors or a combination of them) 
and flows of interest (origin/destination flow volumes, arc flow volumes, route flow volumes, or a com- 
bination of them). Following the classification defined in the recent survey by Gentili and Mirchandani 
[15], two main classes of problems can be identified in this context: (i) locating sensors to fully observe 
flow volumes on the network (Sensor Location Flow-Observability Problems) , and (ii) locating sensors to 
estimate flow volumes on the network (Sensor Location Flow-Estimation Problems). This division derives 
from the observation that the location of sensors on a network (either on the nodes or on the links) can 
be translated into a system of linear equations where the set of variables corresponds to the unknown 
flows and the set of equations comes from the deployed sensors. When the resulting system has a unique 
solution, we say the system is fully observable and therefore, under the assumption of error free data, all 
the flows involved in the system are known (that is they are observable). The resulting location problem 
consists in determining the optimum deployment of sensors on the network that results in an observable 
system (Sensor Location Flow-Observability Problems) and therefore it allows to determine all the flow 
volumes of interest. Observability problems arise, for example, when either path-ID sensors or vehicle-ID 
sensors need to be located on the links of the network to determine route flow volumes or link flow volumes 
(Gentili and Mirchandani [14], [15], [16], Hu et al. [17], He [25], Castillo et al. [2], [3], [4], [5], [6], [7], Minguez 
et al. [22], Ng [23]). On the other hand, when the system is not observable (that is, it is underspecified), 
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it admits an infinite number of solutions. The related location problem consists in determining how to 
optimally deploy sensors on the network so that the derived flow estimates are the best possible. Gener- 
ally, underspecified systems arise when one is interested in determining origin-destination flow volumes 
by locating counting sensors on the links of the network or when there is a budget constraint that limits 
the number of sensors that can be located. This problem has been hugely studied in the literature (see 
for example, Chootinan et al. [8], Eisenman et al. [10] , Elhert et al. [11], Fei and Mahmassani [12], Yang 
et al. [18] Kim et al. [19], Lam and Lo [20], Li and Ouyang [21], Yang and Zhou [26], Zhou and List [27], 
Fei et al. [13], Simonelli et al. [24]). 

In this paper we focus on both the observability and the estimation problems arising when vehicle- 
ID sensors are to be located on the links of the network to determine and/or estimate route flow volumes. 
Both the problems were defined in Castillo et al. [9] and Minguez et al. [22] respectively, where the 
corresponding mathematical formulations were presented and some computational experiments were car- 
ried out to evaluate the quality of the resulting observed/estimated route flows. The experiments were 
performed by solving the proposed mathematical formulations by means of a solver. Since the problems 
are NP-hard only small networks were tested. 

We further study the two problems and provide contributions both from a theoretical and an algorithmic 
point of view. In particular, we improve the existing mathematical formulations by adding some new 
constraints that better define the feasible set of solutions, we also provide efficient solution approaches, 
still missing in the literature, namely two different greedy approaches and a Tabu Search algorithm. The 
efficiency of the proposed approaches is tested on a set of benchmark tests by comparing the provided 
solutions with the optimal ones, when available, returned by the solver CPLEX. 

The sequel of the paper is organized as follows. Next section describes the addressed problems. Sec- 
tion 3 contains the mathematical formulations of the problems and the existing results in the literature. 
Section 4 describes our solution approaches. Computational results are discussed in Section 5. Conclu- 
sions and further research are object of Section 6. 

2 The Vehicle-ID Location Problem: Observability and Estima- 
tion 

In this section we formally describe the addressed problems. We follow the general description given in [15] 
to introduce the Sensor Location Flow-Observability Problem and the Sensor Location Flow-Estimation 
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Table 1: OD Trips, Route Flows, Link Flows related to the network in Figure 1 


Problem when vehicle-ID sensors are to be located on the links of a network. 



Figure 1: An example network with 5 nodes and 8 links. 


Consider the network in Figure 1 with 5 nodes and 8 links, and assume there are 4 different OD pairs 
on the network each connected by one or more routes, namely: the OD pair wi = (1,5) is connected by 
route R\ = {ai, 02 , a.3, 04 } and by route R 2 = {0.1 , 07, a 4}, the OD pair w; 2 = (1, 4) is connected by route 
i?3 = {01, 06) a 8 , 03}, OD pair W 3 = (3,2) is connected by route R 4 = {03,04,05,01} and the OD pair 
W 4 = (4,3) is connected by route R 5 = {04,05,01,02}. Table 1 lists all the OD pairs, the corresponding 
OD trips, all the 5 routes, the corresponding route flow volumes, and all the link flow volumes. For 
example, the OD trips corresponding to OD pair w\ are equal to 27 units, 15 of them use route i?i and 
the remaining vehicles use route R 2 ] flow on link 02 is equal to 37 units and it is equal to the sum of 
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the flow of the routes that use that link, that is routes R\ with flow 15 and route R 5 with flow equal to 
22. The information shown in the table are given here for clarity of the exposition; we assume to know 
only the set of routes, while all the other information, namely, OD trips, link flows, and route flows are 
not known. Our aim is to locate Vehicle-ID sensors on the links of the network to derive the route flow 
volumes (listed in Table 1) that are not known. The location of a vehicle-ID sensor on a link of a network 
allows to univocally recognize a vehicle traveling on the network. When vehicle-ID sensors are located 
on the links we obtain flow volumes of the corresponding links. However, not only flows of the links can 
be obtained, but also, more detailed information can be derived when the same vehicle is detected on 
different links at different times. For example, let us suppose to locate three vehicle-ID sensors on our 
network example, say on links a 4 , a 3 and 014 (we say that links 01, CI3 and 04 are scanned). Since link ai 
is used by all the routes in the network, then the following equation can be derived: 

{link on) f Rl + f R2 + f Rs + f Ri + f R5 = 66 (1) 

where 66 is the intercepted flow on link ai, and f Ri denotes the unknown route flow volume of route R % . 
Similarly, for links CI3 and 0,4 , we can derive the two following equations: 

{link a 3 ) f Rl + f Rs + f Ri = 32 ^ 

{link a 4 ) f Rl + f R2 + f Ri + f R5 = 56 

Moreover, if we consider jointly links ai and 04, the following two additional equations related to the 
number of vehicles that are recognized on both links ai and 0,4 , can also be obtained: 

{links (ai,a 4 )) f Rl + f R2 = 27 
{links (a 4 ,ai)) /r 4 + Jr 5 = 29 

Indeed, when a vehicle is detected, the detection time is also known. Hence, if the same vehicle is detected 
first on link di and subsequently on link a 4 , then we can state that such a vehicle is using any route that 

contains both the links in such an order, in this case the two routes R\ and R.2- While if a vehicle is first 

intercepted on link a 4 and subsequently on link ai then we can say it is using either route R.4 or route R$. 
We refer to this type of equation as the joint recognition equation associated with a subset of vehicle-ID 
sensors. Considering all the possible equations associated with the location set U = {ai,a3,a4}, the 
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following system of linear equations will result: 


(link ai) 
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These equations are those associated with some of the possible ordered subsets of these three sensors. In 
particular, the first three equations are associated each with a single sensor, the subsequent five equations 
are the joint recognition equations associated with the possible subsets of two sensors out of the three 
that are located, and, the remaining two equations are the joint recognition equations associated with 
the three sensors together. Note that, not all the possible ordered subsets are considered, indeed there 
does not exist any route in the network that contains the ordered pair (04,03). The same observation 
can be made for the ordered sets containing all the three sensors: (ai, 04, 03), (03, ai, 04), (04, ai, 03) and 

(04, 03, Oi). 

System ( 4 ) has 5 unknown variables and 5 linear independent equations, hence the system admits a 
unique solution and all route flow volumes fa, i = 1 , . . . , 5 can be univocally determined. We say, in this 
case, that the location set U = {ai, a 3 , a 4 } has associated a fully observable system of linear equations. In 
this context, one could wonder whether there exists another subset of three or less sensors such that the 
associated system of equations is fully observable and hence allows to determine all route flow volumes. 
The following location problem arises: 

Vehicle-ID Sensor Location Flow-Observability Problem: 

What is the minimum number of vehicle-ID sensors to be located on the links of the network and where 
to locate them, such that the associated system of linear equations is fully observable, that is it allows to 
uniquely determine all route flow volumes? 

A different question arises when a limited number of sensors, say k, can be located on the network 
and all the systems associated with all the possible subsets of k sensors do not admit a unique solution. 



Indeed, in this case, how can one decide when a location set is better than another? A possible way is 
(see [15]) to define a location-evaluation criterion g that associates with each location set U a value g{U) 
representing a measure of the quality of the flow estimates obtained by the system of equations associated 
with U. Given such a criterion, we can conclude that a location set U is better than another location set 
U' if the value g(U) associated with U is better than the value g{U') associated with U' . Following the 
approach suggested in [22], for the case where vehicle-ID sensors need to be located, we could associate 
with the location set U an evaluation criterion g{U) that equals the total number of routes that can be 
univocally observed through the system associated with U. In this way, we could conclude that between 
two location sets U and U', U is better than U' if g{U) > g(U')- In our network example, let us consider 
the two locations sets Ui = {01,04} and C/ 2 = {03,05}. The following system is associated with the 
location set U\: 

{link ai) f Rl + /k 2 + f Rs + f Ri + f Rs = 66 

{link o 4 ) f Rl + f R2 + f Ri + f R5 = 56 

{links (04,04)) /fii+/fl2 = 27 

{links (04,04)) /fl 4 + /fl5 = 29 

and the following system is associated with the location set C/ 2 : 

{link o 3 ) f Rl + f Rs + f Ri = 32 

{link a 5 ) f Ri + f Rb = 29 (6) 

{links (a 3 ,a 5 )) f Ri = 7 

By solving system (5) route f? 3 can be univocally determined and g{U\) = 1, while by solving system 
(6) we have 3(1/2) = 2 since routes R4 and R 5 can be univocally determined. Hence, we could conclude 
that, according to the defined evaluation criterion g, f/ 2 is better than U\. Hence, given an evaluation 
function g, the following location problem can be defined: 

Vehicle-ID Sensor Location Flow-Estimation Problem: 

If a given number k of vehicle-ID sensors are to be located on the network, where to locate them so that 
the associated evaluation criterion g is optimized? 


Both the described problems are NP-hard as it will be evident in the next section where the mathe- 
matical formulation is presented and the relationship with the set-covering problem is outlined. 



3 Mathematical formulations 


In this section we provide the mathematical formulation of the two problems described in the previous 
section. To this aim we introduce the concept of scanning map (firstly defined in [9]) associated with a 
location set U. Let us return to our example network. Consider again the location set U\ = {oi, 04} and 
let us associate with each route Ri the set of links C of scanned links of the route (referred to in the 
sequel as scanning set associated with the route). That is: C ^ = (01,04), Cr 2 = (04,04), Cr 3 = (ai), 
Cr 4 = (04, ai) and Cr 3 = (04,01). Let us refer to this collection of sets as the scanning map associated 
with the location {01,04}. Analyzing this scanning map we can conclude: (i) if a vehicle is intercepted 
only on link ai, then it is using route R3: (ii) if a vehicle is intercepted first on link ai and later on 
link 04 then it is using either route R\ or Ra; finally, (iii) if a vehicle is intercepted first on fink 04 and 
later on link ai then it is using either route R4 or R5. That is, this scanning map allows to recognize 
univocally only vehicles of routes R3, while for the remaining routes the same conclusion cannot be de- 
rived. Let us consider to locate one additional vehicle-ID sensor on link 03, the new scanning map is: 
Cr x = (01,03,04), Cr 2 = (01,04), Cr 3 = (01,03), C. r 4 = (03,04,01) and Cr 5 = (04,01). The location of 
a vehicle-ID sensor on link 03 permits to distinguish among those vehicles that use route R\ and those 
vehicles using route R2, and, similarly, we can distinguish those vehicles that use route R4 and those 
vehicles that use route R3. From this example, it is clear that, by using the scanning map associated 
with a location set, it is possible to understand which route, among all, is uniquely identified. Based on 
the scanning map definition we can give an alternative definition of the two problems described in the 
previous section, the mathematical formulation presented in [9] and [22] is based on these definitions. 

Formally, let G = (V, A ) be a network and R be a set of routes defined on the network. Given a lo- 
cation of vehicle-ID sensors on a subset U C A of links, we can define a scanning map SM{U) = { Gp t : 
Cr, = UC\Ri,i = 1, 2, . . . , |i?|}, where each scanning set C'/j, is associated with route R,, and contains the 
scanned ordered list of links of the route. If all the scanning sets in the scanning map are not empty and 
are different, i.e. C 7^ Cr } , for each i ^ j = 1 , 2 , . . . , \ R\, then the scanning map allows to determine 
all the route flow volumes and we can define the alternative version of the observability problem given in 
the previous section: 

Vehicle-ID Sensor Location Flow-Observability Problem: 

What is the minimum number of vehicle-ID sensors to be located on the links of a network and where to 
locate them such that the associated scanning map allows observability of all route flow volumes? 


while, if we define the evaluation criterion g(U) as the number of routes that can be univocally de- 
termined through the scanning map associated with U, we can state the following alternative version of 



the estimation problem: 

Vehicle-ID Sensor Location Flow-Estimation Problem: 

If a given number k of vehicle-ID sensors are to be located on the network, where to locate them so that 
the associated scanning map is such that the total number of univocally determined routes is maximized? 


For a matter of clarity of exposition, in the sequel we will refer to the Vehicle-ID Sensor Location 
Flow-Observability Problem as problem PI, and, to the Vehicle-ID Sensor Location Flow-Estimation 
Problem a problem P2. 

Let us introduce two sets of binary variables: variable x a associated with each link a G A that it is equal 
to 1 if a vehicle-ID sensor is located on link a and it is equal to 0 otherwise; variable y ab associated with 
each pair of links a ^ b £ A that is equal to 1 if a vehicle-ID sensor is located both on link a and link b 
and it is equal to 0 otherwise. The following set of parameters are also defined: 

- p° Ri : is equal to 1 if route Ri contains link a and is equal to 0 otherwise; 

- P R(R .'- is equal to 1 if link a is contained either in R, or in Rj (not in both) and is equal to 0 
otherwise; 

- cr R . R .: is equal to 1 if links a and b are both in route R t and Rj but they appear in a different 
order. 


The mathematical formulation of problem PI, is the following [PI]: 


min5> 
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The objective function (7) is the sum of the total number of the vehicle-ID sensors located on the network 
and is minimized. We refer to constraints (8) as covering constraints as they ensure there is at least one 
vehicle-ID sensor located on each route of the network. Constraints (9) are diversification constraints. 
These constraints ensure the scanning map associated with the location of vehicle-ID sensors allows full 
observability of the route flow volumes. Indeed, the sum on the left of the constraint ensures that for 
each couple of routes Ri and Rj the corresponding associated subsets C Ri and C R . are different either 



because they contain at least one different link or, in case they contain the same set of links, they appear 
in different order. Constraints (10) are logical constraints linking the binary variables x a and y a b- If 
y a b is equal to 1 then the constraint forces both x a and Xb to be equal to 1, otherwise, if y a b = 0 the 
minimization of the objective function drives the value of x a and Xb to be equal to 0. The remaining 
constraints (11) and (12) require the variables to be binary. 

The above mathematical formulation is different from the one provided in [9] where the set of diver- 
sification constraints did not take into account the order the scanned links appear in each set C Rt . Our 
diversification constraints (9)-(10) consider this order and hence the feasible set of solutions of the problem 
is now correctly described. In Section 5 the results of the two mathematical formulations are compared 
on different sets of instances. The results show that our modified version of the diversification constraints 
improves the final solution on those instances where the routes in the network are defined as subsets of 
circular paths like, for example, the two routes R\ and R4 in the network of Figure 1 that are subsets of 
the circular path P = {ai, 03 , 04 , as}. 

Let us consider the additional set of binary variables z Ri associated with each route R,, , such that z Ri is 
equal to 1 if route R, is uniquely identified (that is observed) and it is equal to 0 otherwise. Let us also 
define a weight w Ri associated with route Ri to define its relative importance (see for example [22] for 
possible definitions of these weights). The mathematical formulation of problem P2 is then the following 

[P2] : 


max ^2 w Ri z Ri 
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The objective function (22) maximizes the total weight of the routes that can be univocally determined. 
Constraints (14) and (15) differ from those in the previous model [PI] since they ensure the coverage and 
the diversification of route Ri only if the corresponding variable z Ri is equal to 1. Constraints (16) are 



the same as in the previous model. Constraint (17) ensures at most k vehicle-ID sensors are located on 
the links of the network. Variables are ensured to be binary by constraints (18) - (20). Note that, by 
fixing w Ri = t f Vi?,, £ R, the problem looks for the maximization of the total number of routes that can 
be univocally determined. 

As for the previous formulation [PI], the above formulation of the problem differs from the one pro- 
vided in [22] basically for the characterization of the diversification constraints (15)-(16) that, unlike the 
constraints defined in [22] , take indirectly into account the order the scanned links appear in the scanning 
sets. 

In the sequel of the paper, we consider also two weighted variants of the above described problems 
that arise when an installation cost c a is associated with each link. In particular, a variant of problem 
[PI] requires to determine the minimum cost vehicle-ID sensor location that ensures the observability of 
all the routes in the network. The corresponding model formulation is the following [Plb] 

min5>z a (21) 

s.t. (8)-(12) 

If there is an available maximum budget C max for the location of the sensors, then one could be interested 
in locating vehicle-ID sensors on the links of the network to maximize the total weight of the routes that 
can be univocally determined without exceeding the budget limit. In this case, the following model 
describes this variant [P2b]: 


max ^2 w Ri z Ri 

(22) 

RieR 


S.t. ^2 c a x a < C max 

(23) 

(14) -(16), (18) -(20) 

(24) 


4 Solution approaches 

In this section we present our solution approaches, namely two greedy approaches (Greedyl and Greedy2) 
and a Tabu Search algorithm, to solve problems PI, Plb, P2 and P2b. We describe the approaches for 
problem PI and subsequently we will explain how the approaches have been modified to solve problems 
Plb, P2 and P2b. 
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Table 2 : Set of Covering constraints and Diversification constraints associated with the example in Figure 


4.1 Greedy 1 Approach 

Before describing the greedy criterion used by our Greedyl algorithm, let us introduce first some defini- 
tions and notations. Given a subset of links U C A, we associate with each link a € A two different scores 
according to the number of additional constraints that can be satisfied if a is selected. In particular we 
define the two following criteria: 

• the Covering Criterion: Cov(U, a) is the additional set of routes that can be covered if link a is 
selected (that is, the set of additional covering constraints that can be satisfied if link a is selected). 
The resulting score associated with link a is the size of this set, that is \Cov(U,a)\. 

• the Diversification Criterion: Div(U,a) is the set of additional diversification constraints that 
can be satisfied since either the term ff R . R . x a , or the term cr R . R y a b, becomes equal to 1 in the 
summation of the constraints where x a or y a i, are involved, if a is selected. The resulting score 
associated with link a is the size of this set, that is \Div(U,a)\. 

Let us return to our example in Figure 1 for which the entire set of covering and diversification constraints 
of the mathematical formulation [PI] is given in Table 2 . Let us assume link 0,4 is selected, that is 
U = {0,4}, the corresponding value of the above criteria for each of the remaining links are listed in Table 
3 . The covering constraints that are already satisfied by the selection of link a4 are those related to 
routes: R\,R.2,Ra and P5, while the satisfied diversification constraints are those corresponding to the 
pairs: {Pi, P3}, {P2, P3}, {P3, Ra} and |P3,Ps}. Consider for example link 0,2. Link a 2 has a covering 
score equal to \Cov({a4,},a,2)\ = 0 since there is no additional covering constraint that is satisfied if 
a2 is selected, while the diversification score is equal to 5 ; indeed, the constraints related the pairs 
(Pi, P2}, {Pi, P4}, {P2, P5} and {P4, P5} are satisfied since variable x a2 = 1 and the constraint related 
to the pair {Pi, P 5 } is satisfied since variable y 0204 = 1 . 
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\Cov({di}, a»)| \Div{{ai}, ai)\ 


_ai 1 i_ 

a-i 0 5 

03 1 4 ^ 

05 0 4 

Oq 1 0 

ar 0 3 

a 8 1 0 


Table 3: Scores of each link of the network of Figure 1 for U = {04} . 


According to the two criteria we define an initial ordering of the links by ordering them in decreasing order 
with respect to the Covering criterion for U = 0; ties are broken according to a decreasing order of the 
Diversification criterion for U = 0. We refer to this ranking as Initial Ranking. The Initial Ranking for 
our network example is: {ai, 04, 03, 05, 02,05, 07, a 8 } (see the upper most Table in Figure 2 where the val- 
ues of the Covering criterion and of the Diversification criterion are showed for our example when {7 = 0). 

Our first greedy approach to solve problem PI, iteratively selects a link that maximizes the covering 
criterion; in case of ties the link that maximizes the diversification criterion is chosen; in case of ties 
the link having the best position in the Initial Ranking is selected. The algorithm stops when a feasible 
solution is obtained. A re-optimization step is also applied to the returned solution U. This step tries 
to eliminate redundant links. The details of this algorithm are given in Table 4. In the worst case, the 
running time required by this algorithm is 0(|A||i?| 2 ); indeed in the worst case, the algorithm could, for 
each link, check the entire set of constraints. 

Figure 2 shows the steps of Greedyl when it is applied to solve problem PI on our example network of 
Figure 1. 
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Greedyl - PI 

1. Initialization 

Let U <— 0 be the set of selected links. 

Let A be the set of links in the network. 

2. Greedy step 

2.1 While U is not a feasible solution 

find the link a € A\U with \Cov(U,a)\ maximum; 

in case of ties select the one with | Div(U, a)| maximum; 

in case of ties select the one with the best position in the Initial Ranking 

3. Return U 

4. Re-Optimization Step 

U = {a il ,a i2 ,...,a ih y, 

j h; 

While j > 0 
U <- U\{ aij }; 

if U is not feasible then U U U {o^ } 

3 J ~ 1 

5. Return U 


Table 4: Greedyl Algorithm to solve Problem PI. 
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- Scores of each link of the network of Figure 1 for U = 0 . 

- The resulting initial ranking is: {ai, a4, 03, as, 02, 06, 07, as}. 

- The first selected link is ai: U — {ai}. 


- The second selected link is 04: U = {01,04}. 


- The third selected link is 03 : U = { 01 , 04 , 03 }. 

- The current set U is feasible. 

- The re-optimization step does not eliminate any element from U. 

- The final returned solution is optimal. 


Figure 2: Steps of Greedyl algorithm when it is applied to solve problem PI on the network of Figure 1. 
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Greedy2 - PI 

1. Initialization 

Let U 4 - 0 be the set of selected links. 

Let A be the set of links in the network. 

2. Greedy step 

2.1 While U is not a feasible solution 

Find the link a £ A\U such that W(U,a) is maximum; 
Set U <- {a}; 

3. Return U 


Table 5: Greedy2 Algorithm to solve Problem PI. 

4.2 Greedy2 Approach 

Given a subset of links U C A, our second greedy approach to solve problem PI, iteratively selects a link 
a that maximizes a greedy criterion W(U,a ) that takes into account three main components: (i) length 
of the routes (in terms of number of links) that are covered if link a is selected, (ii) number of routes 
that can be univocally identified if link a is selected, and (iii) a measure of the total number of routes 
that still need to be differentiated if link a is selected. The algorithm stops when a feasible solution is 
obtained. The re-optimization step is applied to the returned solution U. The description of the greedy 
function is given next, while the pseudocode of the algorithm is given in Table 5. Let us better describe 
the function W(U,a) that is the core of the greedy criterion at Step 2.1. Let R be the set of the routes 
on the given network and let U C A be a subset of links. Function W(U,a ) favors the choice of those 
links that, on the one hand, improve most the feasibility of the solution, and, on the other hand, leave a 
set of routes that are easier to be differentiated. These aspects are taken into account by considering the 
sum of three weighted components: 


W(U, a) = Pl h(U, a ) + p 2 f 2 {U, a ) + p 3 f 3 (U, a) 


The first component is computed 


h{U,a)= Y, 

RieCov(U,a) 


Avg 

1^1 


where |i?j| is the total number of links contained in route and Avg is the average number of links of 
each route Rj £ R (for the example in Figure 1 we have Avg = ^ = 3.8). This criterion takes 

into account the total number of routes that can be additionally covered if link a is selected. However, 
unlike the covering criterion, it also takes into account the length (in terms of number of links) of the 
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additional routes that can be covered. Indeed, routes with a fewer number of links are less likely to be 
covered, hence, function /i favors the selection of those links that cover shortest routes. 

The second component f 2 (U,a) is the following: 

f 2 (U,a) = \D(U,a)\. 

This component equals the total number of routes that can be univocally identified by the selection of 
link a. This criterion is slightly different from the diversification criterion introduced previously. Indeed, 
Div(U, a) is the set of diversification constraints that can be satisfied when link a is selected, while 
D(U, a) is the set of routes that are associated with a unique nonempty scanning set in the scanning map 
resulting after the selection of link a. 

Finally, the last component is the following: 

fs(U,a)= Y, V\M(C Ri )\. 

C Ri eSM(UU{a}): 

Ri<£D(U,a) 

RieCov($,a ) 

This component considers the scanning map that results after the selection of link a. In particular, for 
each route that is covered by link a, i.e. Ri £ Cov(<D, a) and that is still not univocally determined, that is 
Ri ^ D(U,a), the associated scanning set Cr± is considered. The total number of routes whose scanning 
set is equal to Cr x , that is \M(CR i )\ = \{Rj £ R. : Rj 7^ R % and Cr. = Cjj 4 }| is computed. Function 
fs computes the sum of the square root of the size of these sets. This criterion privileges the selection of 
those links that, if selected, split the resulting scanning map as much as possible. 

These three functions are properly weighted with weights Pi, i = 1, 2, 3 to define the relative importance 
of the corresponding criterion. In our experiments we set pi = 2p 2 , p 2 = 100 and pz = 1. We chose these 
values to favor criterion fi first, then criterion f 2 , and finally to take into account criterion fz in case of 
ties. 

Figure 3 shows the steps of Greedy2 when it is applied to solve problem PI on the example of Fig- 
ure 1. In the example we set p\ = p 2 = p 2 = 1. Tables 6-8 show the scanning maps associated with each 
link and with sets U = 0, U = {ai}, and U = {«! , ciz} respectively. 

In the worst case, the running time required by this algorithm is 0(|A||i?| 2 ); indeed in the worst case, 
the algorithm could, for each link, check the entire set of constraints. 
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SM({ai}) 

SM({a 2 }) 

SM({a s }) 

SM({a 4 }) 

SM({a s }) 

SM({a 6 }) 

SM({ar }) 

SM({a 8 » 

ill 

(01) 

(a 2 ) 

(« 3 ) 

(a 4 ) 





R2 

(«i) 



(a 4 ) 



(“7) 


Rs 

(ai) 


(« 3 ) 



(“e) 


(as) 

Ri 

(ai) 


(« 3 ) 

(a 4 ) 





R5 

(«i) 

(ct 2 ) 


(a 4 ) 

(a 5 ) 





Table 6 : Scanning maps associated with the singleton sets {a*}, Va, £ A. These scanning maps define 
the scores of each link as shown in the upper table of Figure 3. 



SM({ ai },a 2 ) 

SM({ ai },a 3 ) 

| SM({ ai },« 4 ) | SM({a 1 },a s ) 

I SM({ ai },a 6 ) 

SM({»,},« 7 ) 

| SM({ ai },a 8 ) 

Ri 

(ai,a 2 ) 

(aii 03) 

(ana 4 ) 

(ai) 

(ai) 

(ai) 

(ai) 

R2 

(a 1) 

(01) 


(ai) 

(“1) 

(ai, 07) 

(“1) 


(«i) 

(ai, a 3 ) 

(«i) 

(«i) 

(°ii “el 

(of) 

(ai,a 8 ) 

Ri 

(«i) 

(a 3 , ai) 

(a 4 , ai) 

(a 5 ,ai) 

(ai) 

(ai) 

(ai) 

R n 

(ai, a 2 ) 

(«i> “ 3 ) 

(a 4 ,ai) 

(a 5 ,ai) 

(ai) 

(ai) 

( 51 ) 


Table 7: Scanning maps associated with the sets {{oi},aj}, Va* £ T\{a-| }. These scanning maps define 
the scores of each link as shown in the middle table of Figure 3. 



SM({ai,a 4 },a 2 ) 

SM({ai,o 4 },o 3 ) 

SM({ ai , o 4 },o B ) 

SM({ai,a 4 },a 6 ) 

SM({ai,o 4 },« 7 ) 

SM({ai,a 4 },a 8 ) 

Ri 

(ai, a 2 , a 4 ) 

(ai,a 3 ,o 4 ) 

(a 1.0.1) 

(ai, a 4 ) 

(ai,a 4 ) 

(Oi,a 4 ) 

r 2 







Rs 

(ai) 

(ai, 03) 

(ai) 

(ai, a 6 ) 

(ai) 

(ai, a 8 ) 


(a 4 , ai) 


(a 4 ,a 5 ,ai) 

(a 4 ,ai) 

(a 4 , ai) 

(a 4 ,ai) 


(a 4 , ai, a 2 ) 

(a 4 ,oi) 

(a 4 ,o B ,ai) 

(a 4 ,oi) 

(a 4 , a 1) 

(a 4 ,oi) 


Table 8 : Scanning maps associated with the sets {{ai, 03 }, a,*}, Va* £ A\{ai , a 3 }. These scanning maps 
define the scores of each link as shown in the last table of Figure 3. 
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- Scores of each link of the network of Figure 1 for U = 0 . 

- The first selected link is 01: U = {au}. 


- The second selected link is d4i {di, 04}- 


- The third selected link is d2 (or also <23): U = {01, d4, 02}. 

- The current set U is feasible. 

- The re-optimization step does not eliminate any element from U. 

- The final returned solution is optimal. 


Figure 3: Steps of Greedy2 algorithm when it is applied to solve problem PI on the network of Figure 1. 
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FindNeighbor( U ) 

1. Let U' «- t/; 

2. Randomly select a £U'\ 

3. For h iterations do 

3.1 Randomly select Ri G Cov(%, a); 

3.2 Randomly select, if it exists, a link d such that a £ U fl Rj/, 

3.3 Set U' <- U'\{a}; 

4. Add elements to U' by applying Greedy2 until a feasible solution U is obtained; 

5. Return U. 


Table 9: The procedure FindNeighbor(U). 

4.3 Tabu Search Approach 

In this section we describe the Tabu Search approach we developed to solve problem PI. Given a feasible 
solution U C A for the problem, we defined the neighborhood Nh(U) of this solution as the set of all the 
feasible solutions obtained from U by deleting at most h elements. Formally, the following neighborhood 
is considered: 

k-switch Neigborhood Nh(U ) 

A subset of links U' £ Nh(U) if and only ifU 1 is such that: \U' fl U\ > |t/| — h, that is the two subsets U 
and W share at least \U\ — h links. 

Note that, a neighbor U' € Nh(U) can be obtained from U by randomly deleting h elements and adding 
new elements to restore feasibility; such a procedure, referred to as FindNeighbor(U), is described in 
Table 9. The value of parameter h must be, on the one hand, large enough to ensure to find a neighbor 
U' e N(U)h such that U' ^ U, and, on the other hand, its value should vary according to the size of the 
instance. After a tuning phase we set h = 5 + /3(\U\), where /3(|f7|) is a random value between 0 and 
The tabu list memorizes moves corresponding to feasible solutions. Ideally, the tabu list should store 
the subset of links of each solution that is contained in the tabu fist. Since memorizing such a detailed 
information could be heavy from a memory point of view, we decided to store in the tabu fist an hash key 
H(U) that is univocally associated with solution U . When the algorithm needs to verify whether a given 
solution is in the tabu fist it looks for its hash key in the tabu list. To give an example, let us consider 
U = {a .^ , a l2 . . . . , cii s } to be a feasible solution composed of \U\ links, where a.^ is the identification ID 
of the link. The corresponding hash key H{U) is obtained by applying the following steps: 

• Order the links in U in non decreasing order with respect to their ID and obtain: < a, )2 < 

’ ’ — 3 \U\ 

• 7 ( 0 ) = 1 
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Iteration i-th of Tabu Search 

1. Apply FindNeighbor(U) and obtain U ; 

2. If | O’ | < 1 17* | then U* ^U,T <-T\J 

3. If |0| > j 17* j and H(U) <£ T and |0| < |17**| then U** <- 0; 

4. If after p iterations the incumbent solution U* was never updated then U** <— U * ; 


Table 10: The i-th iteration of our tabu search algorithm 

• 7 (h) = aiUf. + «27 (/i - 1 ) 

The first step is to order the links in the solution with respect to their ID. This step ensures that the 
hash key associated with a given solution does not depend on the order of the links in the solution. The 
hash key associated with U is a value computed by means of a recursive function 7 (-). The value y(/i) is 
associated with the link in position h, h = l,2,...,|t/|, in the ordered set U: such a value is computed 
according to the value of the recursive function y(/i — 1) associated with the link in the previous position 
h— 1 . The hash key associated with U is the value of the recursive function computed for the link in the 
last position. The values aq and (*2 are prime numbers. If the links were identified by prime numbers 
the resulting hash key would be unique. Of course, due to the size of the instance of the problem such 
an identification for the links is not possible. This results in a key that may not be unique. To limit the 
cases of different feasible solutions associated with the same hash key value, one would choose the values 
of ol\ and a 2 as greater as possible. In our experimentation we set aq = 4523 and a 2 = 31. 

Finally, we chose to implement a tabu list with a fixed length equal to 8. The generic iteration step of 
our tabu search algorithm is described in Table 10 where U is the current feasible solution, U* is the 
current incumbent solution, U** is the second incumbent solution and T is the current tabu list. The 
value of parameter p is fixed in our experiments equal to \U\ * 5. The initial feasible solution to start 
the first iteration of our tabu search is obtained by applying Greedy2. Our Tabu Search terminates if 
no improvement in the incumbent solution is obtained after \N\ * 150 iteration, where \N\ is the total 
number of nodes of the input graph. 

4.4 Adaptation to solve problem Plb 

Our solution approaches described to solve problem PI can be easily adapted also to solve problem Plb 
and therefore to take into account the cost c a of installation of a sensor associated with each link. 

The modifications for Greedyl involve the definition of the Initial Ranking and of the greedy step criterion. 
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The Initial ranking is such that the links are ordered increasingly with respect to the link cost; ties are 
broken according first to a decreasing order with respect to the Covering Criterion (computed for U = 0) 
and then according to a decreasing order with respect to the Diversification Criterion (computed for 

17 = 0). 

The greedy Step 2.1 of the Greedyl-Pl (in Table 4) becomes: 

Step 2.1. While U is not a feasible solution 

- find the link a £ A\U with, maximum; 

- in case of ties select the one with Jtf^ ,a ^ maximum; 

- in case of ties select the one with the best position in the Initial Ranking. 

The modifications for Greedy2 involve the greedy criterion W(U,a), that becomes VV ’W“) . Such a mod- 
ification is also considered in the Tabu Search algorithm. Moreover, Step 3 in the i-th iteration in the 
Tabu Search (in Table 10) is also changed to take into account the new objective function. It becomes: 

3. If c(f>) > c(Z7*) and H(U) <£ T and c(U) < c(U**) then U** SM U. 
where c(U) = J2 a eu c « is the sum of the costs of the links in the set U. 

4.5 Adaptation to solve problem P2 and P2b 

The adaptation of the approaches to solve problem P2 and P2b requires to take into account a new 
objective function (the maximization of total weight of the routes that are observable) and the added 
budget constraint involving either the total number of sensors that can be located (problem P2) or the 
total cost of the sensors that can be located (problem P2b). All the changes are described for problem 
P2b and reported below. The same changes can be applied to solve problem P2 by setting C max = k 
and c a = 1, Vo. £ A. Let z(U) define the total weight of the routes that can be univocally determined if 
vehicle-ID sensors are located on the links in U . 

The modifications for Greedyl algorithm involve the Initial Ranking and the Greedy criterion. When 
solving problem P2b the Initial Ranking takes into account the cost associated with each link and it is 
the same we considered for problem Plb, that is, the links are ordered increasingly with respect to the 
link cost; ties are broken according first to a decreasing order with respect to the Covering Criterion 
(computed for U = 0) and then according to a decreasing order with respect to the Diversification 
Criterion (computed for U = 0). The greedy criterion to select link a is modified to take into account 
both the total number of additional constraints that can be satisfied, and the weights of the additional 
covered routes. In particular, the algorithm selects at each iteration a link a such that: (i) the budget 
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Greedyl - P2b 

1. Initialization 

Let U -f- 0 be the set of selected links. 

Let A be the set of links in the network. 

2. Greedy step 

2.1 While ( c(U ) < C max ) and (A\U ^ 0) 

find the link a £ A\U with . 1 _|_ maximum 

in case of ties select the link with the best position in the Initial Ranking; 
if ( c(U ) + c(a)) < C max then A <- ,4\{a} and U <- U U {a}; 
else A <- A\{a} 

3. Return U 


Table 11: Greedyl Algorithm to solve Problem P2b. 


constraint is not violated and (ii) the following quantity is maximum 

\Cov(U,a)\ + \Div(U,a)\ T Ri& Cov(u,a)™Ri 
c(a) w 

where the first element is the ratio between the sum of the covering and diversification scores over the 
cost of the link; while the second element is the ration between the sum of the weights of the additional 
covered routes divided by the average route weight w = — — -■ The pseudocode is given in Table 11. 
Modifications of the Greedy2 algorithm involves the stopping criterion and the values of the parameters 
Pi, i = 1,2,3 associated to the three components that are involved in the greedy criterion. The new 
values assigned to parameters pi are: p\ = 100, p 2 = 2 p\ and ps = 1; this setting favors the choice of 
those links that diversify most. Greedy2 stops when the given budget limit is reached. Finally, the same 
modifications apply to the Tabu Search algorithm. Moreover, Step 3 in the <-tli iteration in the Tabu 
Search is also changed to take into account the new objective function. It becomes: 

3. If z(U) < z(U*) and H{U) £ T and z(U) > z{U**) then U** U. 


5 Experimentation 

In this section we show the experimental results to test the performance of our approaches, in terms of 
solution quality and running time, and to evaluate when our proposed enhanced mathematical formu- 
lations provide a different optimal solutions than the formulations presented by Castillo et al. [9] and 
Minguez et al. [22]. 

Since there do not exist benchmark instances in the literature, we generated two different classes of 
instances. The first one, referred to as Class 1, is a set of instances where the set of routes do not show 
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the circular pattern behavior described in section 3. The second class of instances, referred to as Class 
2, contains instances where the set of routes in the network show a circular behavior. These instances 
are useful to evaluate the effectiveness of our enhanced mathematical formulations in providing better 
solutions with respect to those provided by the existing mathematical formulations. 

Class 1 and Class 2 instances are downloadable from the authors’ website [1] On each instance, we 
compare the results returned by the mathematical models and by our Greedy 1, Greedy 2 and Tabu 
search algorithms for each of the four problems. The greedy algorithms and the Tabu Search algorithm 
were coded in C and run on a 3.4 Ghz Intel i7 processor. The mathematical formulations were coded in 
AMPL and solved by means of the solver CPLEX 11.2. We fixed a threshold on the execution time of 
the CPLEX solver equal to 2 hours. 


5.1 Instances Description 

Class 1 set of instances are modified square grid graphs, where arcs are added randomly between pairs 
of vertices of the graph. In particular, a number of arcs equal to 10% of the total number of arcs of a 
regular grid were randomly added to each grid. We generated square grids whose dimensions are: 8x8, 
10 x 10 and 15 x 15, for a total number of nodes equal to n = 64, 100 and 225, respectively. For each 
grid we generated a fixed number \R\ of routes such that \R\ = n, 2 n, 3 n, 4 n, 5 n. Hence, we generated 15 
different scenarios. Costs on links were generated randomly from a uniform variable and they range from 
1 to 5. Weights on routes were generated randomly from a uniform variable and they range from 10 to 
100. For each scenario, 5 different instances were generated. 

Class 2 set of instances are generated considering networks with a total number of nodes equal to 
n = 40,80,160,240,320,400 and corresponding number of OD pairs equal to j . A fixed number \R\ of 
routes was generated for each OD pair with |f?| = 4,8, 10, 16. We then generated 24 different scenarios. 
Costs on links were generated randomly from a uniform variable and they range from 1 to 5. Weights 
on routes were generated randomly from a uniform variable and they range from 10 to 100. For each 
scenario, ten different instances were generated. 

5.2 Results 

All the results are showed in Tables 12-18. Since computational times and memory requirements to solve 
the mathematical models hugely increase with the size of the instances, we report the results of our mod- 
els for the first 40 instances of Class 1 and the first 60 instances of Class 2, being the time limit reached 
for problem P2 and problem PI for higher dimensional instances. The detailed results on each instance 
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on these limited set, returned by the mathematical models, can be downloaded from the authors’ website 
[1] . Comparisons of the solutions returned by our heuristic approaches with the optimal solution (when 
available), on these limited set of instances, are reported in Table 13 and Table 14 Finally, comparisons of 
the quality of the solutions returned by our heuristic approaches on the entire set of instances are given 
in Tables 15 - 18. 

To compare the effectiveness of the new set of diversification constraints Table 12 shows, for each problem 
and each class of instances, the percentage of instances where our mathematical formulation returned a 
better solution than the one returned by the models in the literature. Hence, for example, the percentage 
of instances of Class 2 where the solution of our mathematical formulation of problem PI returned a 
different (and hence better) solution that the one returned by the existing mathematical formulation is 
equal to 80%, and on the remaining 20% of the instances the two mathematical formulations returned 
the same solution. The new set of diversification constraints reveals to be more effective in diversifying 
routes related to Class 2 instances than those of Class 1 instances. Note that, the budget constraints for 
problem P2b makes the new diversification constraints not effective at all. 

The analysis of the results in the detailed tables downloadable from the website reveal the average 
gap between the optimum solution is 1% on Class 1 instances and 3% on the instances of Class 2. Our 
proposed models require more computational time to be solved than those existing in the literature. 
This was an expected result being the total number of variables and of constraints greater. Additionally, 
comparing the computational times of the instances in Class 1 and Class 2 with the same number of 
routes, we can observe that higher computational times are required to solve instances of Class 1 than 
those of Class 2, revealing that graphs where the routes show a circular pattern are more easy to deal 
with as it was expected. Among the four problems, problem Plb and problem P2b seem easier (in terms 
of required computational time) to be solved than problems PI and P2. In particular, the time limit to 
solve our mathematical formulation of problem PI, Plb, P2 and P2b on Class 1 instances was reached, 
respectively, on 10, 0, 33, 1 instances out of 40. Moreover, our mathematical formulation was not able 
to find a feasible solution within the given time limit for Problem P2 on two instances (namely, instance 
37 and instance 39) of Class 1. While, the time limit to solve our mathematical formulation of the four 
problems was never reached on the 60 considered instances of Class 2 when solving problems PI, Plb, 
and P2b, while it was reached on 8 instances out of 60 when solving problem P2. This can be explained 
with the fact that having different cost associated with the arcs allows, on the one hand, to obtain better 
bounds that allow to cut off part of the branch and bound tree (for problem Plb), while, on the other 
hand, to prune the tree tanks to the budget constraint (for problem P2b). 

Tables 13 and 14 compare the quality of the solutions provided by Greedy 1, Greedy2 and the Tabu 
Search algorithm with the optimum solution (when available) on the 40 instances of Class 1 and the 60 
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PI 

Plb 

P2 

P2b 

Class 1 

20% 

18% 

50% 

0% 

Class 2 

80% 

82% 

95% 

0% 


Table 12: Comparison of the mathematical formulations: for each problem and each class of instances, 
the percentage of instances where our mathematical formulation returned a better solution than the one 
returned by the models in the literature is given. 


instances of Class2 on which we run the mathematical formulations to optimality. In each table, the 
first three columns report the characteristics of each scenario: scenario ID, the number of nodes (n) and 
the number of routes (r). For each algorithm, the tables report the average values taken on the tested 
instances of the relative gap between the returned solution and the optimum solution, when available. 
In particular, this gap is computed as the ratio UB -&p£ mum f or problem PI and Plb, and °P Um ^- LB 
for problem P2 and P2b, where UB and LB are the value of the bounds returned by our algorithms. 
Whenever at least one instance of the scenario has not been solved to optimality by the solver CPLEX 
within the time limit of 7200 secs., the term NA (Not Available) appears in the tables. In such case, of 
course, no gap value is reported for all the tested algorithms. 

Tables 15 - 18 compare the solutions and the computational time of Greedyl, Greedy2 and the Tabu 
Search algorithm on all the instances of Classl and Class2. In each table, the first three columns report 
the characteristics of each scenario: scenario ID, the number of nodes (n) and the number of routes (r). 
For each algorithm, the tables report the average values for each scenario taken on the tested instances 
(5 instances for Classl and 10 Instances for Class 2) of the bound and of the computational time (in 
seconds) 
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The percentage gap from the optimum of the tabu search is always better than those showed by 
Greedyl and Greedy2 as it is expected. This better performance is, however, paid in terms of greater 
computational time as it is evident by the analysis of tables 15 - 18. The percentage gap of the Tabu 
Search ranges between 0% and 2% when solving problems PI and Plb. These gaps are higher when 
solving problem P2 and P2b. In particular, the gaps for problem P2 are not available for many of the 
instances of Class 1 (table 13) while the gaps for problem P2 on instances of Class 2 varies between 4% 
and 6% (table 14). Finally, the tabu search, when solving problem P2b has a gap that ranges from 8% 
to 61% on Class 1 instances and a gap that ranges from 10 % to 15 % on Class 2 instances. The solution 
returned by Greedy2 algorithm is consistently better than that returned by Greedyl. The two algorithms 
are however comparable in terms of computational times which are always less than one second, but for 
Greedyl algorithm on scenarios 20, 23 and 24 on instances of Class 2 for problem P2 and on scenarios 
16,20,22, 23 and 24 on instances of Class 2 for problem P2b (table 18). 
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Table 18: Comparison of the algorithms on instances of Class 2 for problems P2 and P2b. 




6 Conclusions 


We addressed an important problem in the context of traffic control and management related to the opti- 
mum location of vehicle-ID sensors on the links of a network to derive route flow volumes. We considered 
both the full observability version of the problem where one seeks for the minimum number of sensors 
(or minimum cost) such that all the route flow volumes can be derived and the estimation version of the 
problem that arises when there is a limited budget in the location of sensors. Our contribution is both 
theoretical and algorithmic. We improved the existing mathematical formulations of the four addressed 
problems by better describing the feasible region. Being the problems NP-complete, we proposed three 
solution approaches (still missing in the literature) to heuristically solve the problems: two greedy al- 
gorithms and a tabu search metaheuristic. We tested our approaches on two different set of instances 
to better evaluate the quality of solutions provided by our algorithms and by the mathematical models 
when the set of routes show a circular pattern behavior or not. 

The proposed analysis is developed under the assumption (common to other papers existing in the 
literature [9], [22]) that the information regarding OD routes is correct and complete. An interesting 
development of this research is now focused on the analysis of the sensitivity of the solutions with respect 
to the level of knowledge of the OD routes. 
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